Moving faults while unfaulting 3D seismic images

Unfaulting seismic images to correlate seismic reflectors across faults is helpful in seismic interpretation, and is useful for seismic horizon extraction. Methods for unfaulting typically assume that fault geometries need not change during unfaulting. However, for seismic images containing multiple faults and, especially, intersecting faults, this assumption often results in unnecessary distortions in unfaulted images. We developed two methods to compute vector shifts that simultaneously move fault blocks and the faults themselves to obtain an unfaulted image with minimal distortions. For both methods, we use estimated fault positions and slip vectors to construct unfaulting equations for image samples alongside faults, and we construct simple partial di↵erential equations for samples away from faults. We solve these two di↵erent kinds of equations simultaneously to compute unfaulting vector shifts that are continuous everywhere except at faults. We test both methods on a synthetic seismic image containing normal, reverse, and intersecting faults, and we also apply one of the methods to a real 3D seismic image complicated by numerous intersecting faults.